吸 积 盘 随 机 振荡 光度 功率 谱 密 度 中 低频 准 周期 振 沪 现象 的 
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摘要 采用 含 频率 涨 落 噪声 和 指数 形式 关联 随机 力作 用 的 广义 朗 之 万 方程 模型 描述 黑 
洞 吸 积 盘 的 垂 向 振荡 ， 推 导出 吸 积 盘 随机 振荡 光度 功率 谱 密 度 的 解析 表达 式 ， 并 讨论 了 系 
统 参 数 对 功率 密度 谱 中 低频 率 准 周期 振荡 (LFQPOs) 现象 的 影响 。 研 究 结果 发 现 : 选取 合 
适 的 系统 参数 时 ,功率 密度 谱 曲 线 上 出 现 了 一 个 基 频 和 一 个 二 次 谐 频 的 共振 双 峰 LFQPOs, 
基 频 峰 对 应 的 中 心 频率 为 吸 积 盘 振荡 的 特征 频率 ， 随 机 力 关 联 时 间 决 定 了 基 频 峰 的 高 度 和 
度 ， 频 率 噪声 强度 和 烙 滞 阻尼 只 对 二 次 谐 频 峰 产生 影响 。 结 果 说 明了 吸 积 盘 的 随机 振荡 
型 可 作为 LFQPOs 起 源 的 一 种 解释 。 

关键 词 低频 准 周期 振荡 ， 吸 积 盘 ， 频 率 涨 落 噪声 ， 随 机 共振 
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低频 准 周期 振荡 CLFQPOsO 是 存在 于 黑洞 双星 暂 现 源 的 普遍 现象 ,表现 为 X 射线 光 变 
功率 谱 密 度 曲线 上 的 宽 峰 ， 峰 的 中 心 频率 范围 从 几 mHz 到 30Hz03]。 根 据 峰 的 品质 因子 、 
M 噪声 和 相位 延迟 等 因素 ，LFQPOSs 通常 被 分 为 A、B、C 三 种 。 其 中 ，C 型 LFQPOs 是 黑洞 
ST 暂 现 源 中 最 常见 的 类 型 ， 常 在 低 硬 态 和 硬 中 间 态 被 观测 到 ， 频 率 范围 在 几 mHz 到 10 Hz 之 
= 间 , 但 在 高 软 态 也 能 被 观察 到 , 其 频率 可 达到 ~30Hz。C 型 LFQPOs 伴随 着 较 强 的 平台 噪声 ， 
A Jt Hic t LIE ADL ARITA EERO, EK AEEA JL RERE, HE 
主峰 的 中 心 频率 与 振荡 的 基 频 有 关 ， 附 加 峰 是 指 分 别 在 基 频 的 半 频 和 二 次 谐 频 处 出 现 的 副 
谐 频 峰 和 谐 频 峰 P5。 在 少数 情况 下 ， 谐 频 峰 〈 或 副 谐 频 峰 ) 的 值 与 基 频 峰 相当 ， 形 成 一 种 
特殊 的 “cathedral-like”， 即 双 峰 LFQPOs。 如 : 在 XTE J1550-564 和 XTE J1859+226 两 个 黑 
洞 暂 现 源 功 率 谱 的 相关 频率 处 出 现 了 两 个 强 而 窗 的 峰 [67。 
c 虽然 对 LFQPOs 的 观测 和 研究 取得 了 很 大 的 进展 ， 但 是 其 物理 起 源 一 直 是 一 个 有 争议 
© 的 问题 。 其 中 ,一 个 可 能 的 起 源 是 吸 积 盘 的 振荡 。 在 吸 积 盘 中 ,流体 元 受 扰动 的 回复 力作 用 

时 ， 在 水 平方 向 上 作 惯 性 振荡 ， 垂 向 在 艾 道 面 附近 作 简 谐振 荡 ， 且 不 同方 向 上 的 振荡 是 耦合 
的 。 这 些 振荡 按照 回复 力 的 性 质 可 以 分 为 g- 模 (gravity modes)、p- 模 (pressure modes). c- 
模 Ccorrugation modes) 三 种 类 型 。 研 究 发 现 : 这 些 振荡 被 禁闭 在 吸 积 盘 中 一 定 区 域内 传播 ， 
RU Aris “trapped oscillation”&9。 禁 闭 振荡 的 频率 都 与 黑洞 质量 成 反比 ， 可 能 是 黑洞 暂 现 源 中 
观测 到 的 准 周 期 振荡 。 其 中 ，g- 模 和 p- 模 振荡 频率 都 在 高 频 准 周期 振荡 的 频率 范围 之 内 ， 而 
c- 模 是 低频 振荡 ,被 禁闭 在 内 半径 和 内 部 垂 向 共振 半径 之 间 ， 其 振荡 频率 与 LFQPOSs 频率 相 
S40), kbh, Titarchuk 和 Osherovich 提出 了 吸 积 盘 全 局 的 垂 向 振荡 模型 ， 该 模型 假设 整个 
盘 的 垂 向 振荡 可 以 产生 0.1 Hz 左右 的 持续 LFQPOsP!, 

以 上 吸 积 盘 的 振荡 模型 并 不 能 从 观测 的 光谱 上 拟 合 LFQPOs, HMA C 型 LFQPOs 的 性 质 

及 宽带 噪声 的 形状 和 强度 可 以 看 出 ， 准 周期 振荡 行为 应 该 与 含有 附加 噪声 的 随机 系统 有 关 。 
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因此 ，Harko 和 Leung 等 人 提出 用 广义 Langevin 方程 描述 吸 积 盘 的 随机 振荡 ， 并 通过 数值 
计算 ,得 到 振荡 光度 的 功率 谱 密 度 ， 解 释 了 一 些 观 测 结果 山 雪 。Long 等 人 2016 年 利用 此 模 
型 计算 了 三 种 不 同情 况 下 的 光 变 曲线 和 功率 谱 , 并 通过 帘 律 拟 合 功率 谱 来 获得 不 同情 况 下 的 
光谱 指数 ， 发 现 观测 天 体 光 谱 指 数 变 化 是 由 盘 振 荡 在 不 同情 况 之 间 的 转换 过 程 中 引起 的 。 
对 于 这 样 的 随机 系统 ， 在 一 定 条 件 下 ， 能 产生 随机 共振 现象 。 而 随机 共振 是 乘 性 在 随机 动 
力 系统 中 的 非 线性 现象 ， 已 广泛 应 用 于 天 文 领域 的 研究 叫 。 因 此 ， 我 们 团队 讨论 了 吸 积 
盘 随 机 振荡 光度 功率 谱 密度 曲线 中 的 随机 共振 现象 ”"。 并 认为 曲线 中 的 共振 峰 为 黑洞 双 
星系 统 中 准 周期 振荡 现象 的 物理 解释 。 在 这 些 工作 中 , 模型 中 的 随机 项 仅 存 在 于 粘 淆 系数 和 
外 部 作用 力 ， 吸 积 盘 的 振荡 频率 是 固定 不 变 的 ,在 功率 谱 密 度 曲线 中 没有 发 现 谐 频 峰 。 基 于 
此 ， 本 文 提出 用 包含 频率 涨 落 和 随机 力作 用 下 的 广义 Langevin 方程 模型 描述 黑洞 吸 积 盘 的 
垂 疝 振荡 ,计算 吸 积 盘 随机 振荡 光度 的 功率 谱 密 度 , 讨论 出 现 的 双 峰 随机 共振 现象 ， 为 解释 
C 型 LFQPOs 提供 理论 依据 。 
1 吸 积 盘 随机 振荡 的 功率 谱 密 度 
11 吸 积 盘 振荡 模型 方程 

考虑 一 包围 致密 天 体 的 标准 吸 积 薄 盘 , 在 受到 外 界 随 机 力 和 内 部 粘 滞 力 作用 时 ， 盘 中 单 
位 质量 粒子 的 垂 向 振荡 可 用 广义 的 朗 之 万 方程 描述 0 

d?z 


E e£ raiz = f). (1) 


式 中 ，z RI o, 分 别 表示 粒子 垂 向 振荡 的 位 移 和 频率 ，c 是 吸 积 盘 内 部 的 粘 滑 耗 散 作用 系数 ， 

为 常量 。 外 部 作用 力 f(?) 是 指数 形式 关联 的 湛 流 作用 随机 力 ， 满 足以 下 统计 性 质 : 
(f(D)=0,(f(Wf0)= De ^^^, Q) 

式 中 ， (…) 表示 平均 值 。 D 为 常数 ,表示 随机 力作 用 的 强度 ,4 表示 随机 力 关 联 时 间 的 倒数 ， 


即时 间 关 联 率 ， 其 大 小 和 盘 振 荡 频 率 关 系 会 影响 吸 积 盘 振荡 与 清流 作用 的 耦合 03。 在 此 ， 
我 们 设 定 4 小 于 盘 的 振荡 频率 。 


RES 


假设 吸 积 盘 的 垂 向 振荡 频率 受 外 界 因 素 的 影响 存在 涨 落 ， 其 表示 式 可 以 写成 09]; 
o,=0,+n(t), (3) 
频率 涨 落 (0) 是 高 斯 型 白 噪 声 随 机 量 ， 其 统计 特性 为 : 
(O) =0, MONE) = eo ó(t -t?. (4) 


Roh. ecoRRPSSRHEAGRBUAGdE WosHd4EOEXezlmW. HAA e 值 较 大 时 ， 频 率 噪 
声 会 迫使 功率 谱 密度 曲线 上 的 峰 较 宽 ， 而 降低 峰 的 品质 因子 QIE HS cu 为 吸 积 盘 垂 向 振 
荡 的 特征 角 频 率 ， 其 表示 式 为 0; 


1/2 
_ 2 447x10° Hz | 2- y [1-7/0 + Dr | m 
m 5 n ud | 
Xd, m-2MÍM,.M 和 MM 分 别 为 中 心 天 体 和 太阳 的 质量 ，1, = R, /3Rs 是 吸 积 
内 边缘 的 半径 ， As = 2GM /c* 是 史 西 半 径 ; ‘ou = Rus [Ra 和 Pj eee 分 别 表示 
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盘 的 外 边缘 半径 和 调整 半径 ; 面 密度 指数 yY 可 以 取 3/5 或 3/400, 

为 了 计算 的 简便 ， 假 定 B = 6/2 作为 吸 积 盘 内 部 的 粘 滞 阻 尼 系 数 ， 则 (可 以 写成 ; 
pe yen © 
dt dt” 

假设 两 个 随机 量 是 不 相关 联 的 ， 满 足 (f(1)n(1')) =0 。 通 过 对 方程 (6) 作 拉 普 拉 斯 变换 ， 

可 以 求 得 任意 时 刻 粒子 振荡 的 位 移 和 速度 : 


2(t) = (0) + | H E-F) -ns)z(s)]ds, (7) 
v(t) = Q0) + f HG- Fs) -ns)z(s)]ds, (8) 
l dH 

RP, HER :假设 初始 位 移 和 速度 为 540) = zy, v(0) = V6， 则 位 移 和 速度 的 平均 

fA (z(t)) mvt) A: 
(z(t)) = voH (t), (v(t) = voH (t). (9) 

在 Qu > 的 情况 下 , 方程 的 响应 函数 AOA AOA: 
g . e P . 
H (t) = sin ot, H (t) = ——(o cos ot — B sin ot), (10) 
CO [n) 


RF, o-Q4o-f . 
1.2 ”振荡 位 移 和 速度 的 关联 函数 

根据 Mankin 等 人 给 出 的 广义 朗 之 万 方程 中 的 位 移 和 速度 关联 函数 的 计算 方法 PRI， 定义 
位 移 的 方差 go?(1)， 及 位 移 和 速度 的 时 间 关联 函数 (4,7) c (0, 7) X: 


c^ (t)-(z(0)-(QG»]). (11) 
K (t, c)-«[z(t * 7 )-(z(t +T) [z@)-(z(@)) D. (12) 
K, (t, c)- ([v(t + 7)-(v(t +r) ]EmsCO v (0) . (13) 


根据 随机 量 w(t) 和 fF (0) 统计 过 程 的 独立 性 ， 它 们 满足 以 下 统计 关系 : 
(FENG) =S) = 0. ()2(0)) = 0, 
(zz) m) = (t - 6 X2"(t,)) (14) 
则 位 移 和 速度 的 关联 函数 可 以 表示 为 : 
ure Bana eG f »dtt, 


+f HDH (GENAN) a 


(15) 


eor) ^ f Hamas Cra) dhdt, 


+ [HGB fca n@)o? q - 1)dt, 
当 ! 一 >"oo 时 ， 以 上 公式 变 为 : 


(16) 


x (Tt)=D| [ HG) (Ge ^" "anat, +o| HCHG nd. AD 


K, (=D f I, H(t, H(t, )e ^*^ ldt dt, + o I, H(t)H(t*z)dt,, — (18) 


式 中 ，G ”是 浙 近 状态 下 位 移 方差 的 时 间 齐 次 部 分 ,1 > oot, €, (0) 2 67, YU RTI IJ: 
ap. 2D(2f +A) 


TELLE a 19 
di o (4f -€)(@, * 2BÀ * A?) (19) 
13 ”随机 振荡 光度 的 功率 谱 密度 
吸 积 盘 中 单位 质量 粒子 振荡 的 总 能 量 为 : 
E = ele. (20) 


由 于 外 部 随机 力 、 及 内 部 粘 滞 耗 散 和 频率 涨 落 噪声 的 作用 ， 粒 子 的 总 能 量 会 部 分 损失 ， 
损失 的 能 量 即 为 粒子 振荡 的 光度 上 (1) ， 由 方程 (6) 和 (20) 可 得 到 : 
dE 


L(t) =~ = 28v + nav - f (v . (21) 
f£ t — 00 时 ,光度 的 稳 态 平均 自 相 关 函 数 定 义 为 : 
C, (7) = lim(L( 0) -L +T) [LO - (4401) Q2) 


将 方程 (17), (18) 和 (21) 代 入 方程 (22), 38: 
Crr(z) = lim[45 "(0 (t - 2) (8 (0)) - 88" (Ot - 7)8(0)* 
-2p(9"(r- 2) (f()8() -4B(8(t TO) f (09 7) Q3) 
2B f (t+ 7)8(t + 7) 0 (0) + (ft 7)8( OAOD) 
Ot - cy + )e()b( xc] 
AHO) = [ ACASO naso an. ed - [Ha OE) meyer. 
通过 传 利 叶 变换 ， 方 程 (230 转换 成 振荡 光度 的 功率 谱 密度 PSD(v) 


PSD(v) - | NO (24) 
通过 复杂 的 计算 ， 并 忽略 掉 不 含 频率 v 的 常数 项 ， 可 得 到 PSD(v) 的 表示 式 : 
42 4D'4^ (2aa; + a,a,a, — a;a,) " EO, (ao —2a,8° *3a fo — 40,6"), 
2 


Au PEARY 20) 
E (B - AYa, * 2BÀ 4 Az? v?) D’ Aa, " De wo & o0 B(B +a) 
V2Xa, Bola,-2BAy 26 ^ 2m" (48? c Am v?) 
. O(a, * 2A — Az^v?) D'4a,  Deayo(a; «12 A^) 
J2za, 2Bo(a,-28AY —- 2w(a? - 4872?) 
8/2D' Aq? 


十 
Mz (44? - Ag?) v?! a -487A7) 
(25) 


式 中 : 


a,=0, +4’, a, = Bo(12n’v’ - Ao), 


a, -[16z*v'*32(B? — o? )a?v?^ 4160; 48^ -4z?v?), 

a, = Plog - A^), a, =@(a@, +A’), 

a, = PB[A(B^ +A Vv) +248’ v? -8o? z^ v^ +80"), 

a, -16o(B^-z?v^)(z^v! -eQ), a,-(aj*4z?v?y -16o?z?v?. 
2 ”模型 参数 对 功率 谱 密 度 曲线 中 LFQPOs 的 影响 

根据 解析 表达 式 (25)， 可 计算 不 同 参 数 作用 下 盘 振 荡 光 度 的 功率 密度 谱 ， 并 发 现 曲线 
中 出 现 了 单 峰 和 双 峰 的 随机 共振 现象 。 这 些 随 机 共振 峰 即 为 LFQPOs。 下 面 主要 讨论 系统 参 
r(e, 4, B, Vy)Xt LFQPOs 的 影响 。 
2.4. LFQPOs 随 特 征 振荡 频率 v, 的 变化 
在 我 们 的 模型 中 ， 通 过 计算 发 现 随机 力 的 强度 万 不 改变 谱 线 的 轮廓 ， 因 此 在 以 下 的 分 
析 中 取 D=1， 且 吸 积 盘 垂 向 振荡 的 特征 频率 内 = 中 /2 近似 为 C 型 LFQPOs 的 基 频 。 由 方程 


(5) 可 知 ， 吧 依赖 于 吸 积 稳 的 结构 和 中 心 天 体 的 质量 。 选 择 不 同 的 盘 参 数 ， 使 振荡 频率 在 
LFQPOs 的 范围 内 。 在 其 它 参 数 保持 不 变 时 , 图 1 给 出 了 不 同 特征 频率 下 的 功率 谱 密度 曲线 ， 
均 呈 现 出 与 LFQPOs 相 吻 合 的 低频 截断 霸 律 谱 轮廓 。 


PSD[(rms/mean) /Hz] 


TEN! 和 | paia 


0.01 0.1 1 10 
Frequency(Hz) 


图 1 振荡 频率 取 不 同 值 时 ， 振 荡 光 度 的 功率 谱 密 度 曲线 ， 参 数 为 E= 0.38, 8 =0.14=0.01D=1。 


Fig. 1 The PSD of the stochastically oscillating luminosity for different characteristic frequencies and 


£ =0.38, 8 20.,4 20.01, D «1. 
从 图 1 可 以 看 出 ， 当 Yo 取 不 同 值 时 ， 小 于 截断 处 频率 的 低频 部 分 的 功率 变化 很 小 ， 
高 频 部 分 的 谱 斜 率 基 本 保持 不 变 。 然 而 ， 曲 线 中 的 随机 共振 峰 (LFQPOs) 对 w 十 分 敏感 。 
只 有 当 Y 达 到 一 定 的 值 时 ， 才 出 现 一 个 基 频 峰 和 一 个 附加 的 二 次 谐 频 峰 。 双 峰 的 中 心 频率 


分 别 对 应 的 Ye A 2vo 。 随 着 w 的 减 小 ， 两 个 峰 变 低 、 变 宽 ， 且 基 频 峰 变 化 更 快 。 当 vo 减 小 


到 0.07 时 ， 基 频 峰 几乎 消失 , 在 曲线 上 只 留 下 一 个 较 宽 的 二 次 谐 频 峰 。 随 着 Vo 进一步 减 小 ， 
谐 频 峰 也 越 来 越 弱 ， 直 至 其 完全 消失 。 在 yo = 0.02 时 ， 功 率 谱 密度 随 频 率 的 增加 而 单调 减 
小 ， 不 存在 LFQPOs。 这 说 明 吸 积 盘 振荡 的 频率 越 高 ， 由 随机 共振 出 现 的 准 周期 振荡 现象 走 
明显 。 
3.2 LFQPOs 随 e,2 和 /等 参数 的 变化 

由 方程 (21) 可 以 发 现 : 吸 积 盘 的 振荡 光度 主要 是 由 粘 滞 耗 散 、 随 机 力 和 频率 涨 落 噪 
声 共 同 作用 下 产生 的 , 这 三 种 作用 在 模型 中 分 别 由 烙 灌 阻尼 系数 6、 时 间 关 联 率 4 和 噪声 强度 
2 三 个 参数 控制 。 为 了 讨论 各 参数 对 LFQPOs 的 影响 ， 我 们 取 同 文献 [10] 一 致 的 吸 积 盘 几何 
参数 ， 即 ,=1 ， =10" ，7y=3/5， 人 =3。 天 体 的 质量 选取 被 观测 到 频率 在 0.1-10Hz 
范围 内 C 型 LFQPOs 的 黑洞 候选 体 XTE J1550-564 HÆ., HEHE Remillard 等 2006 年 评 
占 的 该 候选 体质 量 范围 8.4M ~10.8M_D， 取 -9.6。 将 这 些 参数 代入 方程 (5)， 计 算得 到 
吸 积 盘 振 荡 的 特征 频率 为 v = 0.289 Hz 。 当 v 固定 不 变 ，s, 4, 8 取 不 同 值 时 ， 振 荡 光 度 的 功 
率 谱 密 度 曲线 如 图 2 所 示 。 曲 线 中 都 出 现 了 期 望 的 单 峰 或 双 峰 的 LFQPOs， 且 各 参数 对 谱 线 
4 影响 各 不 相同 。 
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图 2 4,e 和 p 分 别 取 不 同 值 时 ， 振 荡 光 度 的 功率 谱 密 度 曲 线 ， 参 数 为 (a) B = 0.1, = 0.38; 
(b) B=0.1,4 =0.01;© 4=0.01s = 0.38 EAA IE LFQPOs 的 中 心 频率 。 
Fig.2 The PSDs of the stochastically oscillating luminosity for different values of 4, e and f:(a) f «0.16 20.38; 


(b) 8 20.14 20.01; ©) 4 2 0.01, = 0.38 , vertical lines mark the central frequencies of LFQPOs. 


图 2 清晰 地 表明 ， 当 吸 积 盘 的 振荡 频率 固定 时 ，LFQPOSs 基 频 峰 的 中 心 频 率 不 随 其 它 
参数 发 生变 化 。 这 些 参数 只 影响 LFQPOs 峰 的 高 度 和 宽度 。 如 : 随机 力 时 间 关 联 率 1 对 基 频 
的 峰 影响 较 大 ， 而 对 二 次 谐 频 峰 几 乎 不 产生 影响 〈 如 图 2(a) 所 示 )， 且 随 着 4 的 增 大 ， 曲 线 上 
移 ， 说 明 随 机 力 的 关联 时 间 越 短 ， 谱 密度 越 大 ， 但 是 其 基 频 对 应 的 峰值 却 在 减 小 ， 峰 变 得 更 
宽 更 平 ， 直 至 消失 ， 此 时 LFQPOs 只 包含 一 个 中 心 频率 为 2y 的 单 峰 。 刚 好 相反 ， 频 率 噪声 
强度 e 和 烙 滞 阻尼 系数 8 却 只 对 二 次 谐 频 峰 有 较 大 的 影响 ， 对 基 频 峰 的 影响 确 很 小 (如 图 2(b) 
和 (c) 所 示 )。 当 增 大 e， 或 减 小 p, 功率 密度 谱 曲 线 上 移 ， 二 次 谐 频 峰 变 得 更 加 尖锐 。 而 se 取 较 
小 值 ， 或 8 取 值 较 大 时 ， 二 次 谐 频 峰 都 会 消失 。 以 上 说 明 随 机 力 关联 时 间 越 短 ， 频 率 噪声 强 
ERK, 阻尼 越 小 , 吸 积 盘 振 荡 与 满 流 作用 、 内 部 噪声 和 铬 滞 作 用 更 容易 产生 共振 , LFQPOs 
现象 越 明显 。 
3 总 结 

本 文 考虑 频率 涨 落 噪声 和 指数 形式 关联 随机 力 的 作用 ， 改 进 了 吸 积 盘 垂 向 随机 振荡 的 
Langevin 方 程 模型 ， 发 现在 振荡 光度 功率 谱 密 度 曲 线 上 产生 了 广义 的 单 峰 和 双 峰 随机 共振 ， 
用 来 解释 包含 一 个 基 频 峰 和 二 次 谐 频 峰 的 LFQPOs。LFQPOs 峰 是 由 粘 灌 耗 散 、 频 率 噪声 和 
外 界 满 流 随机 力 共 同 作用 产生 共振 的 结果 ,而 且 基 频 峰 对 应 的 中 心 频率 为 吸 积 盘 振 荡 的 特征 
频率 ww， 随机 力 关 联 时 间 决 定 了 基 频 峰 的 高 度 和 宽度 , 频率 噪声 强度 和 烙 灌 阻尼 对 二 次 谐 频 
峰 产 生 影 响 。 
研究 结果 证 实 了 吸 积 盘 的 随机 振荡 模型 对 于 描述 X 射 线 光 变 是 近似 有 效 的 ， 能 较 好 地 解 
释 LFQPOs 现 象 。 在 模型 中 ，LFQPOs 的 中 心 频 率 仅 由 吸 积 盘 的 振荡 频率 决定 。 可 以 想象 ， 
如 果 通 过 调整 吸 积 盘 的 几何 结构 , 使 其 振荡 的 频率 在 瞬时 获得 足够 大 , 将 会 产生 高 频 的 准 周 
期 振荡 “HFQPOs )， 在 合适 的 动力 学 条 件 下 ， 可 以 通过 共振 解释 观测 到 的 3:2HFQPOs 对 ， 
这 需要 我 们 下 一 步 建 立 更 加 复杂 的 吸 积 盘 随 机 振荡 模型 ， 去 拟 合 观测 到 的 HFQPOs 。 
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Density of the Stochastic Oscillating Luminosity of Accretion Disk 
Chen Peijie Wang Zhiyun 
(School of Physics and Electronic Engineering, Hubei University of Arts and Science, Xiangyang 441053) 


Abstract: In this paper, we improve the generalized Langevin equation model with 
fluctuating frequency and exponential function correlated stochastic force to describe the vertical 
oscillations of the accretion disk, then derive the analytical expression of the power spectrum 
density of stochastic oscillating luminosity, and discuss the influence of system parameters on low 
frequency quasi-periodic oscillations (LFQPOs) in power density spectrum. The results show that, 
for appropriate parameters, the double-peaked LFQPOs with a fundamental peak and a second 
harmonic peak can be found in the PSD curves, and the central frequency of primary peak 
correspond to the characteristic frequency of accretion disk oscillation V; . Moreover, the 
correlation time of stochastic force determines the height and width of the primary peak, and 
frequency noise intensity and viscous coefficient only influence the second harmonic peak. These 
results confirm that the stochastic oscillation model of accretion disk should be an alternative 
interpretation for the origin of LFQPOs. 
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